Transition to Turbulence in 
Coupled Maps on Hierarchical Lattices 
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General hierarchical lattices of coupled maps are considered as dynamical systems. These models 
may describe many processes occurring in heterogeneous media with tree-like structures. The tran- 
sition to turbulence via spatiotemporal intermittency is investigated for these geometries. Critical 
exponents associated to the onset of turbulence are calculated as functions of the parameters of the 
systems. No evidence of nontrivial collective behavior is observed in the global quantity used to 
characterize the spatiotemporal dynamics. 

PACS numbers: 05.45.-a, 02.50.-r 



I. INTRODUCTION. 

Coupled map lattices are spatiotemporal dynamical 
systems comprised of an interacting array of discrete- 
time maps. Coupled map lattices have provided fruitful 
models for the study of many spatiotemporal processes 
in a variety of contexts Qj. In most of these studies, the 
dynamical processes have been assumed to take place on 
uniform Euclidean spaces. However; because of their dis- 
crete spatial nature, coupled map lattice systems seem 
especially appropriate for investigating physical phenom- 
ena occurring in heterogeneous media. In this respect, 
there are some recent models of coupled maps defined on 
fractal lattices 0, || , and also an investigation of coupled 
maps on a Cayley tree || . 

In this article we consider coupled maps defined on gen- 
eralized hierarchical lattices, or multi-branching trees, as 
dynamical systems. Examples of phenomena where hi- 
erarchical structures appear include DLA clusters, capil- 
larity, chemical reactions in porous media [^) , turbulence 
|j|, ecological systems |?J, interstellar cloud complexes 
|p |, etc. Hierarchical structures have also been studied 
in neural nets, because of their exponentially large stor- 
age capacity ||. Although many hierarchical structures 
found in nature have random ramifications, here we study 
the case of simple, deterministic tree-like arrays. This al- 
lows to focus on the changes induced in spatiotemporal 
processes as a result of the hierarchical structure of the 
interactions in the system. 

In Sec. II, coupled map lattice models for the treat- 
ment of generalized hierarchical lattices are presented. 
As an application, in Sec. Ill we study the phenomenon 
of spatiotemporal intermittency as a mechanism for the 
transition to turbulence in coupled maps on hierarchical 
lattices. The model is based on the one introduced earlier 
by Chate and Manneville for regular Euclidean lattices 



in one and two dimensions (T(|, and which captures the 
essential features of the transition to turbulence in ex- 
tended systems. A discussion is given in Sec. IV. 



II. COUPLED MAPS ON HIERARCHICAL 
LATTICES. 

Generalized trees constitute a class of hierarchical lat- 
tices which can be generated for any constant ramifica- 
tion number k. At the initial level (which we call level 0), 
there is one cell which splits into k branches connecting k 
daughters cells, which constitute the level of construction 
1. Each cell then splits into k daughters, producing k 2 
sites at level 2. This construction continues until some 
level n. There are k m cells at the level of construction 
to; to = 0,1 ... ,n. Thus, each cell in the lattice has 
one parent and k daughters, except for the level cell, 
which has no parent, and for boundary cells at level n, 
which have no daughters. The number of cells lying on 
the boundary is k n . 

Each cell in a tree constructed until level n can be 
identified by an index i — 1,2, ... ,N, where N is the 
total number of cells on the tree, or the system size, given 
byiV= (fc n+1 -l)/(fc-l). 

A cell denoted by the index i at a level to is coupled 
only to neighbor cells belonging to adjacent levels on the 
lattice, i.e., to its parent and to its daughters, which 
we indicate by indexes i p and id, respectively. We do not 
consider interactions between cells belonging to the same 
level. Thus, a discrete diffusive coupling can be defined 
according to the following scheme. If i = 0, then the cell 
is at the level of the network and it is coupled to its 
k daughters, which have the indices i = 1, 2 . . . , k. Cells 
with i > are coupled to cells with indices i p and id, 
such as 
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where the function int means integral part, and 
i d = (ixk)+j; j = l,2,...,k; if i< (N-k n ), (2) 
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and id is not denned for boundary cells, which have 
i > (N ~ k n ). The level cell has k neighbors, inter- 
mediate cells have k + 1 neighbors, and boundary cells 
for which i > (N — k n ) possess only one neighbor. As 
an illustration, Fig. (|l|) shows a hierarchical tree with 
ramification number k = 3 and construction level n = 3. 
The lattice size is N = 40. The indices on each cell are 
indicated. In this example, the cell labeled by i = 11 be- 
longs to level n = 2; its parent has index i p = 3, and its 
daughters are labeled by id = 34, 35, and 36, respectively. 




FIG. 1: Hierarchical lattice with ramification k = 3 and level 
of construction n — 3, showing the indices on the cells. 

The state of each site on the lattice can be assigned a 
continuous variable x, which evolves according to a deter- 
ministic rule depending on its own value and the values 
of its connecting neighbors. The equations describing the 
dynamics of a diffusively coupled map hierarchical lattice 
at level of construction n can be written as 

(a) level cell: 

x t+1 (0) = (1 - fc 7 )/(x t (0)) + 7 E / (**(?)) i ( 3 ) 

(b) fori < (N-k n ), 

x t+1 (i)= [l-(fc + l) 7 ]/(z t «)+ 

7/(^(g) + iE J& /feW); {) 

(c) and for boundary cells with i > (N — k n ), 

xt+i{i) = f {x t {i)) + 1 f {xt{i p )) ; (5) 

where Xt(i) gives the state of the cell i at discrete time t; 
i p and id label the N cells on the lattice; 7 is a parame- 
ter measuring the coupling strength between neighboring 
sites, and /(x) is a nonlinear function describing the lo- 
cal dynamics. The above coupled map lattice equations 
can be generalized to include nonuniform coupling and/or 
varying ramifications within the network. By providing 
appropriate local dynamics and couplings, different spa- 
tiotemporal phenomena can be studied on tree- like struc- 
tures. 



III. AN APPLICATION: TRANSITION TO 
TURBULENCE. 

Spatiotemporal intermittency in extended systems 
consists of a sustained regime where coherent and chaotic 



domains coexist and evolve in space and time. The tran- 
sition to turbulence via spatiotemporal intermittency has 
been studied in coupled map lattices whose spatial sup- 
ports are Euclidean jnj [ll], [l2| and also in fractal lattices 
H . A local map possessing the minimal requirements for 
observing spatiotemporal intermittency is fllOfl 
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if x G [0, 1] 
if x > 1, 



(6) 



with r > 2. This map is chaotic for f(x) in [0, 1]. How- 
ever, for /(x) > 1 the iteration is locked on a fixed point. 
The local state can thus be seen as a continuum of sta- 
ble "laminar" fixed points (x > 1) adjacent to a chaotic 
repeller or "turbulent" state (x £ [0,1]). 

In regular dimensions, the turbulent state can propa- 
gate through the lattice in time for a large enough cou- 
pling, producing sustained regimes of spatiotemporal in- 
termittency [|10|. Here, we investigate the phenomenon 
of transition to turbulence in hierarchical lattices using 
the local map / (Eq.(||)) in the coupled system described 
by Eqs. (3)-(5). The local parameter is fixed at the value 
r = 3 in all the calculations. As observed for regular lat- 
tices, starting from random initial conditions and after 
some transient regime, our systems settle in a stationary 
statistical behavior. The transition to turbulence can be 
characterized through the average value of the instanta- 
neous fraction of turbulent sites F t , a quantity that serves 
as the order parameter |10| , We have calculated (F) as 
a function of the coupling parameter 7 for several hierar- 
chical lattices from a time average of the instantaneous 
turbulent fraction F t , as 



(7) 
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About 10 4 iterations were discarded before taking the 
time average in Eq. (W) and T was typically taken at 
the value 10 4 . Fig. (0) shows (F) vs. 7 for hierarchi- 
cal lattices with ramification numbers k = 2,3 and 4, 
and levels of construction n — 12,8 and 7, respectively. 
We have verified that increasing the averaging time T or 
the lattice size N does not have appreciable effects on 
our results. The standard deviations about each value 
of (F) are very small on the scale of Fig. (|2|). In fact, 
the fluctuations of F t become smaller for increasing ram- 
ification k and/or for higher levels of construction. In 
comparison, high dimensional Euclidean lattices ]l3|, a s 
well as fractal lattices || and globally coupled maps [Q , 
exhibit large oscillations in the time series of the instan- 
taneous turbulent fraction F t for some parameter values. 
In those cases, the fluctuations about the mean values 
of the turbulent fraction do not decrease with increasing 
lattice size or averaging time, pointing to the presence of 
nontrivial collective behavior in the systems. In contrast, 
normal statistical properties are observed in hierarchical 
arrays of coupled maps, even in the cases when these lat- 
tices have the same numbers of local connections than the 
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other arrays. Note that fractal lattices are also spatially 
nonuniform. However, the behavior of a statistical vari- 
able such as the turbulent fraction, is different for hierar- 
chical lattices. Calculations of the mean field in Cayley 
trees using local logistic maps neither show evidence of 
nontrivial collective behavior Qj. Thus, the topology of 
the underlying geometry seems to have a relevant influ- 
ence on the kind of behavior that may emerge in global 
quantities in spatiotemporal systems. 



fc, 0.-1 




FIG. 2: Mean turbulent fraction {F) as a function of the 
coupling 7 for different hierarchical lattices. Thick left curve, 
k = 2, n = 12; middle curve, k — 3, n — 8; thin right curve, 
k — 4, n — 7. The local parameter is fixed at r = 3. 



There exists a critical value of the coupling 7 C for the 
onset of turbulence which becomes smaller for increasing 
ramification of the trees. The threshold for intermittency 
is achieved when the probability of reinjection into the 
unit interval is balanced by the probability of escape from 
this interval, which is constant for fixed r. The former 
is proportional to 7 and to the number of coupled neigh- 
bors, one of them at least being turbulent. Thus, j c must 
decrease with increasing k, for fixed r. Fig. (||) shows 7 C 
as a function of the ramification number k, verifying this 
assumption. 
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As for other lattices, the transition from a laminar 
regime to turbulence can be characterized by critical ex- 
ponents like (3 which scales the variation of the order 
parameter near the transition point: (F) ~ (7 — 7 C )' 3 . 
We numerically checked this relation. Our results are 
given in Fig. (Q) for lattices with different ramification 
numbers, using a log- log plot. Table (|) gives the critical 




FIG. 4: Log-log plot of (F) vs (7- 
k = 2 (+); k = 3 (x); k = 4 (o). 



7c) for hierarchical lattices: 



exponents (3 calculated from the slopes of Fig. (Q) for dif- 
ferent ramifications k, with their respective errors. The 
values of (3 for these hierarchical lattices with different 
ramifications are of the order of that of a one-dimensional 
lattice, which is = 0.250 ± 0.005 (To). Thus, the sta- 
tistical properties of the transition to turbulence for hi- 
erarchical lattices are closer to that of a one-dimensional 
lattice in this scenario. It should also be noticed that the 
global dynamics of one-dimensional coupled map lattices 
neither display nontrivial collective behavior jl5| . 



k 


n 


P 


2 


13 


0.21 ±0.01 


3 


8 


0.22 ±0.02 


4 


7 


0.25 ±0.03 


5 


7 


0.28 ±0.02 



FIG. 3: Critical coupling 7 C as a function of the ramification 
k. 



TABLE I: Critical exponents f3 for hierarchical lattices with 
different ramifications k. The corresponding levels of con- 
struction n are also indicated. 



IV. DISCUSSION. 

We have introduced coupled map lattice systems with 
hierarchical couplings. These models are especially 
suited for investigating spatiotemporal phenomena oc- 
curring in heterogeneous media with an underlying tree- 
like geometry, and they may also represent the hierarchi- 
cal interactions taking place in some organized systems. 

Although we employed the simple Chate-Manneville 
map for illustrating the transition to turbulence on hier- 
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archical lattices, more sophisticated coupled map models 
may prove useful in the study of fully developed turbu- 
lence in fluids, where dissipation of energy occurs in hi- 
erarchical cascades from larger to smaller structures || . 

In hierarchical lattices, as defined in Sec. II, there ex- 
ists a unique path connecting any two given elements on 
the network. This same topological property occurs in a 
one-dimensional lattice with open ends. The similarity in 
the scaling behavior at the transition to turbulence be- 
tween hierarchical lattices of different ramifications and 
a one-dimensional lattice suggests that topological prop- 
erties of the underlying lattice are relevant in the global 
behavior of spatiotemporal systems. 

In addition, the global quantity Eq. (0) follows a nor- 
mal statistical behavior in hierarchical lattices, while 
nontrivial collective oscillations arise in other geometries. 
Hierarchical lattices have different topological properties 
than those of regular higher dimensional Euclidean ar- 
rays or fractal lattices which have previously been con- 
sidered as coupled map systems. No closed loops exist 



in hierarchical lattices, and the number of local connec- 
tions is not uniform since an appreciable fraction of cells 
possessing only one neighbor lie on the boundary. As it 
has been suggested j| [k| [hJ , the topology of the lattice 
seems to play a more important role on the emergence 
of collective behaviors than the number of local connec- 
tions or the dimensionality of the space. The failure to 
observe collective behavior in our models also points to 
this conclusion. Future work on hierarchical and other 
nonuniform lattices may contribute to answer this and 
other remaining questions on the collective dynamics of 
chaotic extended systems. 
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